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We revisit the problem of calculating the probability current for discontinuous states, such that 
may arise in atom trapping or as a result of projective measurements. In the first passage time 
representation, the problem reduces to evaluation of a localised wave originating from the discon- 
tinuity, whose interference with the initial state determines the transfer of probability. Depending 
on the type of discontinuity, the current behaves as t 1 / 2 , t 3 ^ 2 or const(t). Our approach generalises 
earlier work on this subject. 
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I. INTRODUCTION 

Quantum states whose wavefunctions experience, at 
certain points, discontinuity in the value or the first 
derivative may arise in such applications as atom trap- 
ping by means of strong laser beams pQ, or quantum 
measurements aimed at detecting a particle in a given 
region of space (2]. One may wish to estimate the rate 
at which the probability is transferred across a discon- 
tinuity, for example, because it determines the presence 
and the type of the Zeno effect |2]-jB]. This rate can- 
not, however, be obtained by means of the conventional 
current operator owing to the singular behaviour of the 
kinetic energy. In Ref.[2] the authors evaluated the cur- 
rent at the edge of a truncated state by integrating the 
Feynman propagator over the region of support of the 
original wavefunction. The purpose of this paper is to 
generalise their results to the case of an arbitrary discon- 
tinuity which may occur in the initial state. For this we 
will employ the quantum first passage time (FPT) repre- 
sentation [TJ-JTU], leading to the physically appealing pic- 
ture in which a localised wave is emitted from the point 
where the the wavefunction or its first derivative suffers 
a jump, and whose interference with the rest of the state 
determines the probability current at the discontinuity. 
The rest of the paper is organised as follows: in Sect .II we 
introduce the FPT and in Sect. Ill we define the emit- 
ted wave, S~9(x,t), whose purpose is to repair the initial 
discontinuity at finite times. In Section IV we use the 
conservation of probability to derive some useful prop- 
erties of 5'9(x,t). Section V gives the short time limit 
of the probability current and specialises to a number of 
particular cases. Section V contains our conclusions. 



II. EVOLUTION OF THE WAVEFUNCTION IN 
THE FIRST PASSAGE TIME REPRESENTATION 

For a particle of mass M in one dimension, we wish 
to evaluate the probability current for a complex valued 
wavefunction ty(x) = ^i(x) + i^2(x), at a point which 



we choose to be the origin, x = 0. The particle moves in 
a potential V(x) so that the Hamiltonian takes the form 
(we use h = 1) 

H = -d 2 x /2M + V(x), (1) 

and without loss of generality we choose the ^(x) to van- 
ish at some a < and b > . The task is not triv- 
ial since we will assume that both the wavefunction and 
its first derivative may be discontinuous at x = 0, i.e., 
(/(±0) = lime -> 0/(±e)) 

*(-0) + *(+0), *'(-0) ^ *'(+0), (2) 

where the prime denotes differentiation with respect to 
x. Thus, we have 

■*(x) = * L (x)9 L (x) + * R (x)0 R (x), (3) 

where 9 L (x) = 1 for x < and otherwise, and 
9 R {x) = 1 — 6 L (x). To analyse the development of, say, 
\^ L ) = J dx^ L (x)\x) we invoke the first crossing time 
decomposition of a general evolution operator discussed 
in H-HD] 

exp(-i£ft)|$) = Pcxp(-iPHPt)P\<f>) + (4) 

-i [ dt 1 exp[~iH(t~t 1 )][H,P}Pexp(-iPHPt 1 )P\^) 
Jo 

which describes the time evolution, with a Hamiltonian 
H, of a state initially localised within a subspace H-p 
onto which projects the projector P, P\$>) = |$). Equa- 
tion Q has the standard interpretation [TU]: the system 
moves within %-p until a time t\ when it leaves the sub- 
space for the first time, after which it may or may not 
reenter it. The first term in the r.h.s. of Eq. Q corre- 
sponds to the scenario in which the system has not left 
H-p by the time t. 

The choice of the projector P appropriate for our pur- 
pose requires some care. We begin with the real part of 
V L (x), %i(x), define a Y so that *f(0) + ai*f(0) = 0, 
and then choose P ai = J2 n l^nX^nl witn H\4>n) = 
En\4>n)i subject to the boundary conditions </>^(a) = 



2 



0, 0^(0) + cti<Pn (0) = 0) an d the additional require- 
ment that <fin( x ) = 0, for < a; < 6. Since the eigenfunc- 
tions of a real Sturm-Liouville problem form a complete 
orthonormal basis, the operator P ai is indeed the pro- 
jector on a functional space containing our state of inter- 
est, tyf(x). With this, the reduced evolution operator in 
Eq.dll takes the form 



where erfc(z) — (2/y/n) f°° exp(— z 2 )dz is the com- 
plimentary error function |13] . It is readily seen that 
5^(x,t), continuous with discontinuous first derivative, 



it 



2Mtt' 



5*'(+0) = -^'(-0) = 1/2, (9) 



P ai exp(-iP ai Hp ai t)P ai = 



(5) 



serves to repair initial discontinuities in the value and the 
first derivative of 1 $ l (x)6l(x) at x = 0, as illustrated in 
Fig.l. Indeed, from Eqs.Q and (|8| for t > we have 



The commutator [H, P ai ] is independent of ai, as 
shown in the Appendix A, namely (f\[H,P ai ]\g) — 
-[/*' (0)5(0) - /*(0).g'(0)]/2M which provides the last 
ingredient required in Eq.Q. Repeating the above steps 
for imaginary parts of the left wavefunction, ^2 (x), in- 
troduces the parameter a 2 , *a (°) + "2*2 '(0) = 0. We 
then have an exact result 



(x\exp(-iHt)\^ L ) = tpi(x,t\ai) + iip 2 (x,t\a 2 ) + ^ x 



dt x {d x ,G'(x, x' = 0,t - ti)[y»i(0, t\a x ) + i<p 2 (0, t\a 2 )} 



-G(x, 0, t - hM(0, t\ai) + iip 2 (0, t\a 2 )}} (6) 

where (fi, 2 (x,t\ai t2 ) = (x\U ai 2 (t)\^>i 2 ) are the results of 
evolving the real and imaginary parts of the wavefunc- 
tion on an interval [a, 0] with the Hcrmitian Hamiltonians 
corresponding to (possibly) non-physical boundary con- 
ditions at x = 0, <fi, 2 (x 7 t\ai. 2 ) = — ai^tpi 2 (x, t\a\ :2 ). 
and G(x,x',t) = (x\ exp(— iHt)\x'), 




III. SMOOTHING OF THE DISCONTINUITY 

Evaluation of the current at t = only requires 
the short time limit of the Eq.(|6|, which allows us 
to put U ai (t) « U a2 (t) R3 1, in the integral and 
replace the terms in the square brackets by <3/(— 0) 
and 0), respectively. Further, for a small t, 

the full propagator may be replaced by the free- 
particle one [TT],[T2] G(x,x',t) « G (x - x',t) = 
y/M/2mt exp[iM(x-x') 2 /2t] , with which Eq.^ reduces 
to [V L (x,t) = (x\exp(-iHt)\W L )] 



FIG. 1. (colour online) Short time evolution of the ground 
state of an infinite well of a length xo (thick solid): a) 
truncated at three quarters of its length (thick dashed), for 
t/t = 0.001, to = Mxq; b) with the right wall of the well re- 
moved (thick dashed), for t/to = 0.05. Also shown are various 
terms in Eq.|7|. 

f L (+0,t) = * L (-0,i) = *(-0)/2 + <J*(0,t)*'(-0) +0(t) 



<J L (+0,t) = * L (-0,t) = *'(-0)/2 + 5#"(0)*(-0) +0(t) 

We note also that 6^(0, t) is highly oscillatory and 
rapidly decays away from the origin, 

3/2 
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'2 , nix 2 



exp{iMx^/2t). (10) 



V L (x,t) pa 6 L (x)[$ L (x) - i(x\H\* L )t] 
+^(-0)S^'(x, t) + &(-0)6y(x, t). 



(7) 



We note that S^(x, t) is superposition of all waves emit- 
ted by a constant point source at the origin for all 

< ti < t, 



Finally, changing x — ¥ —x, we can repeat the calculation 
for Adding the two contributions then yields 

V(x, t) w *(x) - i[6 L (x)(x\H\V L ) 
-\8^'{x, t)A* + 5V(x, i)A*'], 



■9 R (x){x\H\9 R )]t 
(11) 



6V(x,t) = - 



2M 



G (x,ti)dti 



it 



— exp(iMx 2 /2t) + ^-erfc{^Mx 2 /2it), (8) 



A* = *(+0) - *(-0), A*' = #'(+0) - *'(-0) 

where the second term is responsible for the Schroedinger 
current in a continuous state (see below), while the last 
term describes an additional wave emitted at the point 
of discontinuity. 
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IV. CONSERVATION OF PROBABILITY 

The evolution of the wavefunction is unitary, which 
implies jl \^{x,t)\ 2 dx = J° \t> L {x)\ 2 dx + J* \V L (x)\ 2 dx. 
With a discontinuity characterised by ^(±0) and ^'(±0), 
there are (discounting an overall phase) seven indepen- 
dent real parameters in Eq.(lll. Inserting Eq.([T7|) into 



the above yields, therefore, seven conditions involving in- 
tegrals of S^ix, t) and S^'(x,t), two of which we will 
consider next. Since 5 1 $(x, t) is highly oscillatory except 
in the vicinity of x = [cf. Eq.( |10[ )], we may expand 
$! L > R (x) w *(t0) + *'(±0)a; and extend the integration 
to ±oo. Choosing ^ R {x) = and #'(-0) = Eq.|TT]) 
then yields 



/>oo />0 

/ \5V'(x,t)\ 2 dx = -Re / 5V' 

Jo J — oo 



xl/2 

(x)dx = (12) 



1-JM-k 



Similarly, for ty R (x) = and ty(x) vanishing at the ori- 
gin, \&(— 0) = 0, we have 



/>oo />oo 

/ \6^(x,t)\ 2 dx = Re x6^(x)dx 
Jo Jo 



,3/2 



6VM 3 tt 



(13) 



which can, if one wishes, be verified directly (see Ap- 
pendix C). In the next Section we will also need the re- 
lations 



5^{x,t)dx = - / x5V'(x,t)dx = , (14) 

o Jo ^Mi 



which cannot be deduced from conservation of probabil- 
ity, but are readily obtained by integrating Eq.([8|. 



V. THE SCHROEDINGER CURRENT 

The probabilty current at the discontinuity, x — 0, is 
defined as 



J(t)=0 t \^(x,t)\ 2 dx = d t P H (t) (15) 
Jo 

with the probability P R (t) determined by the interfer- 
ence between various terms in Eq.(lll. Retaining the 
principal contributions, we have 



P*(t)« f \* R (x)\ : 
Jo 



dx 



{[|*(-0)| 2 - |*(+0)]| 2 7T'e / 6V'(x,t)dx 



2|#(-O)||*(+O)|sinA0iIm / 5V'(x,t)dx} (16) 

Jo 

poo 

2i?e[(**'(+0)A*- **(+0)A*') / 6^(x,t)dx+ 

Jo 



/•OO 

A**' A* / 5f*(x,t)5t>'{x,t)]dx} 
Jo 

±{[l*'(-0)| 2 -|*'(+0)| 2 ]i?e / xm(x,t)dx 



-2|*'(-O)p'(±O)|sinA0 2 Jm / xSVx, t)dx} 

Jo 



arg[tf(+0)] - arg[*(-0)], A0 2 = 



where A(f>i 

arg[^ / (+0)] -arg[^ r/ (-0)]. Above we have used Eqs.([l2l 
( 13 1 and (fl4| , which indicate that in Eq. fl6|) the terms 



in the first, second and third curly brackets are propor- 
tional to i 1 / 2 , t and i 3 / 2 , respectively [13]. Thus, to the 
leading order in i,we have: 

A. Continuous state. If no discontinuity is present, 
A*' = A* = 0, Eq.([l6| reduces to the first term in 
the second curly bracket and, as t — > 0, the current J{t) 
takes the usual form, 



J(t) 



2M 



[** (0)tf (0) - ** (0)W'(0)] ~ const(t) (17) 



B. Discontinuity in the value of the function. For 
A* ^ 0, regardless of the values of A*', *(+0) and 
\l/'(±0), Eq.(16l is dominated by the terms in the first 
curly bracket, and the current behaves as t~ x l 2 (see also 
Ref.0), 



J(i)«{[|*(-0)| 2 -|*(+0)] 
-2|tf(-0)||tf(+0)|sinA<£i} 



(18) 



Ay/Mnt 



-1/2 



Note that if the modulus of the initial wavefunction is 
continuous, the current is proportional to the sine of the 
jump in the phase, Afii. 

C. Discontinuity in the first derivative, and ^(0) 7^ 0. 
For A* = 0, A*' ^ and V R (x) ^ the leading con- 
tribution to the current comes from the first two terms 



in the second curly bracket of Eq.(16), and using Eq.(14l 
we have 



J(t) « -Im{**(0)[*'(+0) + #'(-0)]/2} - const(t), (19) 



which reduces to Eq.([17] if *'(+0) = W(-0). 
D. Discontinuity in the first derivative, and ^ (0) = 0. 
For A* = 0, A*' ^ and * fl (0) = from Eqs.JlT) and 
( 26 1 we have 



J(0«{[|*'(-0)| 2 -|*'(+0)| 2 ]/4 + 
|*'(-O)p'(+O)|sinA0 2 /2} 



(20) 



TrAf 3 



,1/2 



Note that if the right hand space is empty, ^> R (x) = 0, 
Eq.(20) reduces to the result of Ref.0. 



VI. CONCLUSIONS AND DISCUSSION 

In summary, a discontinuity present in a wavefunc- 
tion if?(x) at t = 0, can be seen, for t > 0, as a source 
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of an additional wave propagating on both sides of the 
point where Vt'(x) and/or ty'(x) experience jumps. At 
short times, the wave consists of the amplitude emitted 
by a constant point source, S^(x,t) in Eq, ([Sj, and its 
spacial derivative, weighed by the jumps in the values 
of *i>'(x) and the ^(x), respectively. These additional 
terms serve to repair initial discontinuities, and their in- 
terference with the rest of the wavefunction determines 
both direction and rate of the probability transfer. If 
the value of wavefunction is discontinuous, the probabil- 
ity is always transferred at a rate ~ y/t, and the current 
behaves as t^ 1 / 2 . In the special case where the modu- 
lus of iff(x) is continuous, the current is proportional the 
the sine of the jump in its phase. As suggested in (2] 
such states can be obtained as a result of non-detection 
of atom, which eliminates, partially or completely, its 
wavefunction in a specified region of space. For a discon- 
tinuity in the first derivative at a point where ty(x) does 
not vanish we recover the constant current given by the 
usual expression (17 1 with ^'(x) replaced by the mean 
[tf'(-0)+tf'(+0)]/2[cf. Eq,Q]. Such states would arise 
after an instantaneous collapse of a narrow (5-) barrier 
of a finite transparency. Finally, for ^(0) = 0, e.g., in 
a state obtained by switching off a laser induced atomic 
trap [l],[6], the current increases as i 1 / 2 and contains a 
term proportional to the sine of the jump in the phase of 
<&'(z). 

This work was supported by the Basque Government 
grant IT472 and MICINN (Ministerio de Ciencia e Inno- 
vaciUn) grant FIS2009-12773-C02-01. 



and 



[4> n {x)6 L {x)}g"[x)dx 



(22) 



= / g(x)<f%(x)dx +[g'(0)+g(0)/a}(t> n (0). 



With -<P»(x)/2M + V(x)<t> n {x) = E%<f> n (x) and 
^ n <p n (0)(j) n (x) = 5(x), we find that the terms 
J2n {f\<t>n)E n {(j) n g) and f*(0)g(Q)/a cancel and the com- 
mutator (f\[H,P a ]\g) = -[f*' (0)0(0) ~ fW(0)]/2M 
is independent of the choice of a. 

VIII. APPENDIX B. AN ALTERNATIVE 



DERIVATION OF EQ.(ll) 



Writing ^(x,t) = ty(x) + $>(x,t), and solving the 
Schroedinger equation for $ using the Green's function 
technique yields 

/•oo 

®(x,t) =*(x) -i / dt'0(t-t')x (23) 



dx'(x\ exp[-iH(t - t')}\x')(x'\H\^} 
Differentiating Eq.(|3| we have 
(x\H\^ R ) = e L , R (x)[-^ L ' R "(x) + V{x)^ R {x)] 

±^ I [26(x)^ R '(x)+S / (x)^ R (x)} (24) 



Inserting, Eqs.(24l into (23) and sending t — > we have 
Eq.dlil. 



VII. APPENDIX A. COMMUTATOR OF P a 
WITH THE HAMILTONIAN. 



We need to evaluate matrix elements of the 
commutator (f\[H,P a ]\g) = E„{</ l-%»X<Kls> " 
{f\4> n ) {4>n\H\g)} Differentiation and integration by parts 
yields (assuming f(a) = g(a) = 0), 



r(x)[4> n {x)6 L {x)]"dx 



(21) 



f*(x)^(x)dx + [/*'(0) + /*(0)/c# n (0) 



IX. APPENDIX C. VERIFICATION OF EQ.(Ti3|). 



Inserting Eq.fl8| into ( 13 1 and performing the Gaussian 



integration over x yields 

/•OO 

h{t) = I \SW(x,t)\ 2 dx = (25) 
Jo 

(1/8VttM 3 ) f dh ( 1 dt 2 (h - t 2 y 1/2 = t 3/2 /6V^M3 
Jo Jo 

On the other hand, xGo(x,t) = (tfiM)G' (x,t) so that 



h{t) = / x6V{x,t)dx = (1/2M 2 ) / dt't'G (0,t') (26) 
Jo Jo 



t 3 / 2 /3V2mAP 



and I\ = Rel 2 - 
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